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Abstract. We study the non-equilibrium time evolution of the average transverse magnetisation and end- 
to-end correlation functions of the random Ising quantum chain. Starting with fully magnetised states, 
either in the x or z direction, we compute numerically the average quantities. They show similar behaviour 
to the homogeneous chain, that is an algebraic decay in time toward a stationary state. During the time 
evolution, the spatial correlations, measured from one end to the other of the chain, are building up 
and finally at long time they reach a size-dependent constant depending on the distance from criticality. 
Analytical arguments are given which support the numerical results. 

PACS. 75.40.Gb Dynamic properties (dynamic susceptibility, spin waves, spin diffusion, dynamic scaling, 
etc.) 

05.70.Ln Non-equilibrium and irreversible thermodynamics 
05.30-d Quantum statistical mechanics 



1 Introduction 

Non-equilibrium properties of quantum spin chains at very 
low or even zero temperature have been the subject of 
several investigations through decades. Some precursor 
studies were done by Niemejer [Q and Tjion || at the 
end of the sixties. They have shown that the relaxation 
is algebraic instead of exponential as predicted by Ter- 
wiel and Mazur using a weak coupling limit || . More re- 
cently, several authors have considered the relaxation of 
inhomogeneously magnetised initial states and obtained 
analytical results for various classes of free fcrmionic spin 
chains f§,[|,|]|7]|| ■ It was pointed out in ref. |jj that at some 
special wave vectors, characterizing the inhomogeneity of 
the initial state, a slowing down of the relaxation can occur 
for free fermionic models that have a gapped excitation 
spectrum. In the context of aging phenomena, two-time 
non-equilibrium functions were considered in ref. ||pj for 
the Ising and XX quantum chains. Aging occurs for in- 
termediate times in such systems. All these out of equi- 
librium studies were done for systems without any disor- 
der although alternating and dimerized models were also 
considered previously Q]. In this work we study the in- 
fluence of random spin couplings on the non-equilibrium 
relaxation properties. For that purpose we consider the 
paradigmatic random Ising quantum chain. 

The paper is organized as follows: in section 2 we 
present the random Ising quantum chain and discuss some 
of its main equilibrium properties. After the canonical di- 
agonalization, we express the time-dependent operators 
we are interested in. In the next section, we first give the 
general expressions for the non-equilibrium quantities we 



are considering, namely transverse magnetisation and end- 
to-end correlation functions for two different initial states 
and give the numerical results and analytical analysis of 
this quantities. The last section summarizes and discusses 
our results. 



2 Random Ising quantum chain 

The Hamiltonian of the random Ising quantum chain is 
given by 



(1) 



where the as are the Pauli matrices and </& and are ran- 
dom independently distributed couplings and transverse 
fields. The system possesses a quantum phase transition 
at zero temperature controlled by the parameter 



5 = [Inh] - [In J] 



(2) 



where [.] means an average over the disorder distribution. 
The quantum control parameter d separates the param- 
agnetic phase (5 > 0) from the ferromagnetic one (5 < 0). 
Exactly at S — the system is at the critical point. This 
system has been studied intensively by numerical and an- 
alytical means. In particular, using a RG procedure ini- 
tialy introduced by Ma, Dasgupta and Hu jlO) , Fisher has 
obtained asymptotically exact results (confirmed by nu- 
merical works) on the static properties and equal time 
correlations in the vicinity of the critical point |pT[. In 
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addition to Fisher's results, new numerical and analytical 
results have been obtained by various methods and map- 
pings away from the critical point 12, l^ , p^ , p^[ , in the so 
called Griffiths phase [|l6Ul7j] . In particular, the behaviour 
of the singular quantities in the Griffiths phase are all 
characterized by a single dynamic exponent z(S), which is 
continuously varying with the control parameter S |l4| , |l5t . 
More recently, Igloi has extended the exact RG treatment 
of Fisher into the Griffiths phase on both sides of the crit- 
ical point 

The Hamiltonian (|l|), with free boundary conditions 
and L sites, is readily diagonalized after a Jordan- Wigner 
mapping M 



n-l 



i=l 
n-1 



(3) 



(=i 



and a canonical transformation in terms of Fermi opera- 
tors r\ and r/ + pjj: 



(4) 



where the excitation energies e q are the positive solutions 
of the 2L x 2L eigenvalue problem 



TV q = e q V q 



(5) 



where T is the tridiagonal matrix 



/ hi 
h x Ji 
Ji h 2 



h 2 



V 



'■ Jl-i 
Jl-i h L 
h L / 



(6) 



The 2L components eigenvectors V q are usually split into 
two L components eigenvectors <j) q and ip qi such that 
V q (2k - 1) = -<p q (k) and V q {2k) = %{k). 

Stationary and time-dependent expectation values of 
spin operators are expressed easily in terms of the Fermi 
operators r\ q , T] q . Working in direct space, it is more conve- 
nient to expand these quantities in terms of the operators 
A and B related to the diagonal fcrmions through 



q 

B n = ^2^J q (n) (77+ ■ 



(7) 



These operators satisfy the anticommuting algebra 
{^4n ; ^4m} = 2<5„ irri , {B n: B m } = — 25 n ^ m and {A n: B m } 



0. From the time evolution of the diagonal Fermi opera- 
tors, ri+(t) = e im r\^e~ lHt = 77+ e 16 "*, one obtains the time 
evolution of the ^4s and Bs: 

An(t) = ^2(A n A m ) t A m + (A n B m ) t B m (8) 

m 

and 

B n (t) — ^^{B n A m ) t A m + (B n B m ) t B m (9) 

m 

where the basic contractions (..)t, proportional to the 
time-dependent anticommutators, are given by 



(A n A m ) t = ^2 cos(e q t) (10) 

q 

(B n B m ) t = ^2 ^qi^qi 771 ) cos(e g t) (If) 
q 

and 

(A n B m ) t = (B m A n ) t = i ^ <t>q( n )'<Pq( m ) sin(eqt) . (12) 

q 

For the homogeneous chain at the critical point, these con- 
tractions are simply expressed in terms of integer Bessel 
functions J9|,||. In the disordered situation, we evaluate 
them numerically. 

Utilising the previous expressions, one finds for the 
time evolution of the of = B1A1 Pauli operator at site I 

a f (*) = {( B iBk)t(AiA k ) t - (BiA h ) t (AiB k ) t )ai 

k 

+ J2({BiB t ) t (A l A 3 ) t - (BiAj)t{AiBi} t ) BiAj 

+ Y / (BiA l )t(A l A J ) t A l A J +Y / (BiB l ) t (A l B J ) t B l B J . 

(13) 

Contrary to the operators of, the operators of the form 
BiA m are string operators in terms of the Pauli matrices. 
For example we have for the terms BiA m : 



B l A m > l+l = (-l) L + m+1 af ( J] -I I < ( 14 ) 

so that in a general state it is hard to evaluate the ex- 
pectations (<P\BiA m \i&) . In the same way, the end-to-end 
correlation operator erf (^crf^), which gives insight of the 
development of the correlations in the chain during the 
time evolution, is given by 



with 



°i(*K(t) =A 1 (t)B L {t)Q 



2 = n (-°$ 
3=1 



(15) 



(16) 



the charge operator which is a conserved quantity since 
[H, Q] = 0. Inserting the developments M) and (||) into 
the expression (|l^) , one obtains the final form of the end- 
to-end correlation operator. 
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3 Off equilibrium dynamics 
3.1 Expectation values 

In what follows, we are interested in the non-equilibrium 
relaxation properties of the random Ising quantum chain. 
The time evolution of the system is governed by the 
Schrodinger equation and is given formally by the unitary 
evolution 

\&(t)) = exp(-iHt)\<P) (17) 

since energy is conserved. The initial states we consider 
in this work are eigenstates of the Pauli matrices, that is 
either |sg) or \z) such that crf' z \x, z) = \x,z). These states 
are easily accessible experimentally by applying a strong 
magnetic field in the desired direction. The choice of such 
initial states is not only motivated by their physical ac- 
cessibility but as well because they permit close analytical 
formulae for the magnetisation profile and end-to-end cor- 
relations Consider first the \z) state. In this state, 
the a x - y operators act as flip operators and from formu- 
lae (H) and (Q the non- vanishing expectations are only 
those of the a z terms, such that finally we obtain 

m' s {l,t) = {z\o' l (t)\z) = 

{{BiB^tiAiA^t - {BiA k )t(AiB k ) t ) ■ (18) 

k 

Utilising Q\z) = (— l) L |z), the end-to-end correlation 
function in the z state for a chain of even site number 
is given by 

C' L {t) = (z\*f(t)of(t)\z) = 

Y,{{MB k ) t {B L A k ) t - {AiA k ) t {B L B k ) t ) . (19) 

k 

If the initial state is the |x) state, the role of the flip op- 
erators are played by the a v,z operators. Then, the con- 
tributing terms into ( [l3| ) are only the BiAi + \ = afaf +1 , 
whose expectation value is 1 in |x). This leads to 

ml(l,t) = (x\of(t)\x) = 

5^ ((B,B fc ) t <^i4 fc+1 ) t - {BiA k+1 ) t {AiB k ) t ) . (20) 
k 

The connected end-to-end correlation function C£(i) = 
(x\af(t)al(t)\x) - (x\af(t)\x)(x\al(t)\x) is simply 

Cf / (t) = \(A 1 B L ) t \ 2 . (21) 

The simplicity of the formulae obtained so far is due to 
the locality in terms of the A and B operators of the non- 
vanishing expectations entering here. For other quantities, 
for example expectation of af at a general site I, the for- 
mula is no more local and one should use Wick's theorem 
in order to evaluate the string operators appearing. This 
leads to huge determinants that are not efficiently calcula- 
ble numerically, especially if one has to perform averages 
over disorder realizations. In this case, a direct space cal- 
culation is probably more appropriate. 



3.2 Results 

In the numerical study, we use the following distribution 
for the couplings and fields: 

<J) = {\ J i [M (22) 
v ' ID otherwise 

and 

p(h) = S(h - h ) ■ (23) 

The critical field h c is obtained from the vanishing of the 
control parameter S, that is h c — e _1 . Average quantities 

1 N 

[^ = A7^ Q(l) (24) 

i=l 

are computed using typically N ~ 50000 disorder realiza- 
tions for chains of size up to L ~ 100. 

In figure 1 we show the time dependence of the 
average surface and bulk magnetisation, [m^(l,t)] and 
[m z z (L/2,t)} respectively, at the critical point starting 
with the \z) state. For long times, the non-equilibrium 
surface and bulk transverse magnetisations reach a size- 
independent constant value m s z ,b (oo). The relaxation to- 
ward m s z ' b (oo) is clearly slower than exponential and pre- 
sumably algebraic, with superimposed oscillations with 
period of about 2n. This new time-scale is related 
to the disorder distribution properties. We have tested 
the robustness of those oscillations against disorder dis- 
tributions. For exponential like distributions, P{J) = 
aexjp(-aJ), they have a period r varying linearly with 
the width parameter a, meaning that the microscopic time 
scale r is set by the average coupling [ J] , as r ~ [ J] ~ 1 = a. 
So one expects the presence of these oscillations for every 
disorder distributions with a finite first momentum. 




0.2 1 1 1 1 1 1 1 

10 20 30 

t 



Fig. 1. Surface and bulk non-equilibrium average transverse 
magnetisation at the critical field h c starting with the \z) state. 
The thick line corresponds to the fit discussed in the text. 
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In the case of the surface magnetisation, the relaxation 
behaviour is very well fitted by the form i [cos(1.2t — 
1.46) — 0.43 cos(0.3t — 1)] + 0.5 as shown in figure 1, which 
valid the power law relaxation behaviour. For the bulk 
magnetisation such a fit is much more hazardous since 
the oscillations have a more complicated structure. 
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Fig. 3. Asymptotic values of the non-equilibrium average 
transverse magnetisation as functions of the distance to the 
critical field h /h c . Lines with symbols refer to the disordered 
situation while thick lines to the exact pure case solutions. The 
solid(dashed) line corresponds to the |a;)(|z)). 
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Fig. 2. Surface and bulk non-equilibrium average transverse 
magnetisation at the critical field h c starting with the \x) state. 



Figure 2 shows the same quantity starting with the 
\x) state at time t = 0. The surface and bulk non- 
equilibrium transverse magnetisations approache alge- 
braically the constants m*' & (oo). The features observed 
for the average magnetisation are very similar to the ho- 
mogeneous critical case where one finds m z z ' x (t) = 1/2 ± 
Ji(4t)/4t, with the + (-) sign for the z(x)-state @g§. 
Thus the relaxation toward 1/2, is characterized by the 
power law i~ 3 / 2 . 

The asymptotic values of the non-equilibrium mag- 
netisation are depending on the value of the field h Q . We 
present the dependence of these quantities on the trans- 
verse field in figure 3 for increasing sizes for the disor- 
dered case. We give also, in order to compare, the exact 
asymptotic values for the pure case |]||. The behaviour 
of the asymptotic averaged magnetisation doesn't seem to 
have any singularity near the critical field contrary to the 
pure case. Nevertheless, finite-size effects should not be 
excluded too quickly in the vicinity of the critical point 
also the curve shown in figure 3 have almost collapsed. 
One may notice that as in the homogeneous situation, the 
asymptotic limits take the same value at the critical field 
for both initial states. 

We turn now to the discussion of the average non- 
equilibrium end-to-end correlation, [C^ z (t)]. Figure 4 
show the behaviour of this quantity for various fields. 
Since the velocity of the excitations are finite, the end- 
to-end correlations are vanishing up to a threshold time 
t(L) cx L. No signal starting from one end has yet reached 
the other end. For t > t(L), the correlations start to build 
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Fig. 4. Average end-to-end correlation functions at the critical 
field h c and L = 20 starting with the \x) state. The dashed line 
is given by (^). Insert: same as before for the \z) state. 



up and finally for very long times they reach a finite value 
C^' z (oo) depending on L. In figure 5 we present the size 
dependence of the asymptotic value C£(oo) for h D = h c , 
h > h c and h a < h c . 

The asymptotic behaviour of the end-to-end correla- 
tors can be understood as follows. To fix the ideas, con- 
sider the simple connected end-to-end correlator C£(i) = 
\{AiBi,)t\ 2 ■ It can be split into a time-dependent quantity 
and a time independent part: 



Cl(t) = F(L,t)+G(L) 



(25) 



where G(L) = 1/2 
T,q^ q > r qq <[cos(e q - e q ,)t - 



- cos(e q + e, 



and F(L,t) 
i)t] with r qq > 
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Fig. 5. Size dependence of the asymptotic value [G{L)\ for 
h — h c (circles), h — 0.8h c (squares). Inset : Size dependence 
of the asymptotic value C£(oo) for h = 1.2h c (triangles). The 
solid lines indicate the analytical predictions. 



(f) q (l)(f) q i(l)ip q (L)ip q r(L). For t> 1, the arguments of the 
cosines are random phases and the F(L,t) contribution 
vanishes. In fact, since the bottom spectrum is behaving 
at the critical point as e q ~ exp(— const. L 1 ^ 2 ) p3|, the 
phases will be random for t > e _1 ~ exp(const.L 1 i 2 ). The 
average correlator is thus asymptotically given by 



[Cl(t » 1)] ~ [G(L)} 



(26) 



In the sum over all the modes, the dominant contribution 
comes from the lowest mode which is localized near the 
surface, the others being exponentially smaller. Thus, for 
large sizes one has 

[G{L)\ = \ « . (27) 



Physically, for a first excitation t\ vanishing faster than 
1/L, 4>i (1) and ip\{L) give the equilibrium magnetisation 
in the x direction at both boundaries [|22|]23|1. One has 



(of) = 0i(l) = 




-1/2 



(28) 



and a similar expression for (erf) with hj/Jj replaced by 
hh-j I Jh-j- At the critical point, for a typical sample the 
fluctuations of the couplings generated along the chain 
lead to 

Eir(^) 2 - ex p( cii/2 ) (») 

i j=l \ j / 

thanks to the central limit theorem, where c is a positive 
constant. So that the contribution to the average [G(£)] 
of a typical sample is of order exp(— cL 1 / 2 ). Stronger con- 
tributions to the average [G(L)] come from rare samples 



where the sum J^f 1 T[j=i (hj/Jj) 2 is of order one. As no- 
ticed in several previous works [ p^p^ , |25[ | we can relate this 
problem to the surviving probability of a one-dimensional 
walker making L steps with an absorbing boundary at the 
origin. The quantity [G(L)] is related to the probability 
Pl(0 < yi < ul) that a walker starting near the origin, 
on the positive side, ends after L steps at a position y^ 
larger than all previous positions yi without visiting neg- 
ative sites. This is due to the fact that in [G(L)j it is 
the product of both right and left magnetisations that en- 
ter. For 6 = 0, that is at the critical point, the walker 
is unbiased and from the walker interpretation one has 
[G(L)] ~ 1/L. On the contrary for 6 > there is a drift 
proportional to 6 toward the absorbing boundary. In this 
case, contributions of order one are very rare and the av- 
erage [G(L)\, proportional to the surviving probability, 
scales as L~ z l 2 exp(-aL) with a > a constant This 
form fits very well the datas as it could be seen in the inset 
of figure 5. For 6 < 0, that is in the ferromagnetic phase, 
the walker is drifted off the absorbing boundary and for L 
sufficiently large, the typical fluctuations of order L 1 / 2 are 
not sufficient to make the walker cross the origin. [G(L)] 
is reaching a finite constant for L large enough. It is inter- 
esting to notice again here that it is the connected part of 
the average end-to-end non-equilibrium correlator [C£(f)] 
that reaches a constant for t ^> 1. These three different 
regimes are summarized in figure 5. Same arguments can 
be applied to explain the similar behaviour of the average 
correlator [C£(t)]. 



4 Summary 

We have studied numerically the non-equilibrium relax- 
ation properties of the random Ising quantum chain by 
calculating the average transverse magnetisation and end- 
to-end correlation functions. Starting with a fully magne- 
tised initial state, the random system relaxes toward a 
stationary state. The relaxation of the average transverse 
magnetisation, toward a constant value depending on the 
transverse field, is clearly not exponential due to the pres- 
ence of too large oscillations, but most probably algebraic. 
For the same reasons it is hard to extract from the numer- 
ics the exact algebraic decay. This behaviour looks like the 
decay of the pure model without any disorder. Neverthe- 
less, contrary to the pure case where quantum interference 
plays a central role, in the disordered system the averaging 
process over many samples is predominant. 

The time-evolution of the end-to-end correlations, 
which gives some simple insight of the development of 
the spatial correlations in the quantum chain, is some- 
how more instructive. For times smaller than a threshold 
time t(L), they vanish since both ends are not causally 
connected. After this limiting time, they start to grow 
up and for very long times they finally reach a constant 
value depending on the transverse field and on the sys- 
tem size. For the \x) initial state, utilising a random walk 
argument, one can explain the scaling of this stationary 
size-dependent value. At the critical field, the main contri- 
bution to the average correlation comes from rare samples 
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that are strongly correlated due to large domains of strong 25. D. Karevski, R. Juhasz, L. Turban and F. Igloi, Phys. Rev. 
bonds, giving a contribution of order one, while the contri- B 60, 4195 (1999) 
bution of a typical sample is exponentially small with the 
system size. So the average value is completely determined 
by the distribution of rare samples. This leads to a l/L 
scaling behaviour of the stationary value. In the ordered 
phase, the contribution of strongly correlated domains is 
very important and gives rise to a size-independent sta- 
tionary value, as confirmed by the drifted walker anal- 
ogy. On contrary, for a large transverse field, the walker 
is driven toward the absorbing boundary, as explained in 
the previous section, and the stationary value behaves as 
exp(— cL). In conclusion, contrary to the pure chain where 
after a sudden jump to a value L~ a at t = r(L), the 
threshold time |9|, the end-to-end correlation decreases 
toward zero with a stretched exponential like behaviour, 
in the disordered case the asymptotic behaviour is sim- 
pler. This discrepancy lies again in the fact that quantum 
interferences in the pure case are important while they 
play no significant role in the averaging process in the dis- 
ordered situation. The stationary behaviour is just related 
to the statistical distribution of large strongly correlated 
domains. 
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